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A theoretical calculation is performed for the ultrafast spin dynamics in nickel using an exact 
diagonalization method. The present theory mainly focuses on a situation where the intrinsic charge 
and spin dynamics is probed by the nonlinear (magneto-)optical responses on the femtosecond time 
scale, i.e. optical second harmonic generation (SHG) and the nonlinear magneto-optical Kerr effect 
(NOLIMOKE). It is found that the ultrafast charge and spin dynamics are observable on the time 
scale of 10 fs. The charge dynamics proceeds ahead of the spin dynamics, which indicates the 
existence of a spin memory time. The fast decay results from the loss of coherence in the initial 
excited state. Both the material specific and experimental parameters affect the dynamics. We 
find that the increase of exchange interaction mainly accelerates the spin dynamics rather than the 
charge dynamics. A reduction of the hopping integrals, such as present at interfaces, slows down 
the spin dynamics significantly. Besides, it is found that a spectrally broad excitation yields the 
intrinsic speed limit of the charge (SHG) and spin dynamics (NOLIMOKE) while a narrower width 
prolongs the dynamics. This magnetic interface dynamics then should become accessible to state of 
art time resolved nonlinear-optical experiments. 

PACS numbers: 78.47.+p, 78.20.Ls, 75.70.-i 



I. INTRODUCTION 

Recently ultrafast spin dynamics in ferromagnetic met- 
als has attracted a great deal of attention due to its 
possible applications, such as the ultrafast magnetic 
gates. The experimental observation was that, upon 
the excitation of a femtosecond laser pulse, a sharp 
decrease of the magnetization occurs on a time scale 
of 100 fs 0-§, which is far beyond the characteristic 
time scale of spin-lattice interaction. Similar results 
have been independently found in pump and probe lin- 
ear magneto-optics, nonlinear magneto-optics and two- 
photon-photocmission. However the interpretation of 
this behavior has not been given on the same footing 
and remains somewhat speculative. To reconcile these 
intriguing results from different experimental processes 
such as linear and nonlinear optics is rather crucial and 
is one of the goals of the present theoretical study. More- 
over the demagnetization is very similar to the conven- 
tional one H. 

In the SHG experiment J2|, the M(T) curve is estab- 
lished after the electron thermalization is finished while 
the electron and lattice have not reached a common equi- 
librium yet, which indicates a purely electronic feature 
of the ultrafast spin dynamics. Traditionally the spin- 
lattice coupling sets the speed limit of the demagneti- 
zation process, typically on 100 ps ||, but here this is 
clearly not the case. Our previous theoretical studies 
clearly demonstrated that the dephasing of initial states 
is the origin of the spin dynamics on the femtosecond 
time scale [[j]^]] , as probed by transient reflectivities and 
linear magneto-optics. We found that the intrinsic speed 
limit is about 10 fs. This mechanism is a pure quantum 



effect, resulting from the interplay between band struc- 
ture and electron correlation. In this paper, we focus on 
the nonlinear optical response of spin dynamics. We take 
a Ni monolayer as an example. 

This paper is arranged as follows. In Section II, we 
discuss our theoretical scheme while the main results are 
given in Section III. Finally, the conclusions are given in 
Section IV. 



II. THEORETICAL SCHEME 

It is well-established that in ferromagnetic transition 
metals, in particular in Ni, the electron correlation plays 
an important role even in the ground state and possesses 
a significant impact on the excited states, as evidenced 
e.g. by the famous photoemission satellite structure [||. 
This becomes especially true in the nonlinear optical 
process on the ultrafast time scale where highly excited 
states are frequently involved. Therefore, we employ an 
exact-diagonalization method which explicitly avoids a 
perturbative treatment of electron correlation. Within 
our scheme, one does not need to introduce any damping 
term to obtain a correct dephasing time. Our previous 
results showed that the typical time scale is basically set 
by the dispersion of the bands and the strength of elec- 
tron correlation. We begin with a generic Hamiltonian 
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where Us. 



is the on-site electron interaction, 



which can be described in full generality by the three pa- 
rameters Coulomb repulsion U, exchange interaction J, 
and the exchange anisotropy A J H. The set of param- 
eters used for Ni is given in Ref. ^0|. c\ a (c$ CT ) are the 
usual creation (annihilation) operators in the orbital i 
with spin a (a =|, j). £v{K) is the single-particle energy 
spectrum for band v of the nickel monolayer. n va (K) is 
the particle number operator in momentum space. Hso 
is the spin orbit coupling. Since this is a typical many- 
body problem, one cannot solve it without simplifica- 
tion. In order to obtain a tractable model, we first build 
a two-hole basis set. Within this basis set, for each Ni 
atom, the dimension of the Hilbcrt space is 66. The ma- 
trix elements of electron correlation for each atom can 
be obtained analytically. For each K point, the elec- 
tron correlation is embedded in the crystal field as given 
by the bandstructure. This treatment of correlations is 
analogous to a frequency dependent self-energy correc- 
tion. Within this simplification, we are able to exactly 
diagonalize the Hamiltonian for each K point explicitly 
In order to characterize the spin and charge dynam- 
ics clearly, we calculate both these intrinsic quantities: 
S z (t) = <*(0)|&|*(t)) and N(t) = (¥(0)|#|¥(t)>, and 

(2) 

the nonlinear (magneto-)optical susceptibilities Xxzz and 
Xzzz- Here S z — \{n-\ —raj.), N = (n T + rij), which are 
directly related to the observable NOLIMOKE and SHG 

(2) 

yields, respectively. Since \xxzz((^^t)\ and \S z (t)\ mainly 

(2) 

reflect the spin response while \ X lVz(u,t)\ and |JV(i)| re- 
flect the charge response, they will be used as indicators 
to evaluate spin and charge evolutions, respectively. We 
find: 
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where \kl > is the eigenstate with the eigenvalue E k r, 
p(E kh t) =< ^(t)\kl >. 



III. RESULTS 

Before we come to our main results, we would like to 
demonstrate that our Hamiltonian can reasonably de- 
scribe some basic experimental results. It has been well- 
established that a prerequisite to acquire a ferromagnetic 



ground state is a nonzero Coulomb interaction U and ex- 
change interaction J. We can simply check this by setting 
both U and J to zero. Doing so, we find that the ground 
state is a singlet, i. e. paramagnetic state, which contra- 
dicts the ferromagnetic nature of nickel. This proves the 
importance of U and J. Once we use the generic sets of 
U and J of nickel, we obtain a triplet as its ground state, 
from which we can calculate the magnetic moment, 0.88 
]iB- This magnetic moment is larger than that in the 
bulk material, which is consistent with the experimen- 
tal observation. In this case, also the satellite structure 
well-known from photoemission experiments appears in 
the spectrum quite naturally. 



In Fig. 1, we firstly show the effect of exchange cou- 
pling J on the \x X zz(v,t)\ and \XzM^,t)\ as a func- 
tion of time t. The probe frequency lu is fixed at 2 
eV. The initial state is prepared to be 2 eV above the 
ground state with a Gaussian broadening as large as 20 
[ eV, which opens almost all the possible decay chan- 
nels. Such a large distribution width corresponds to a 
very large laser spectral width. As previously explained, 
this choice aims at revealing the real intrinsic speed limit 
of the spin dynamics in our system, which is then not de- 
layed by experimental constraints. In Figs. 1(a) and 
1(b), the generic parameters of Ni monolayers are used. 
There are several interesting features that should be men- 

(2) 

tioned. One notices that in Fig. 1(a) |xa;zz(w,t)| first 
comes up very quickly and reaches its maximum at about 



2-3 fs. Then \ X 



(2) 



, t) | undergoes a sharply decreas- 



ing envelope and oscillates with a very short period. 

(2) 

The dynamics of \xxzz(oJ,t)\ settles down around 10 fs 
(decay to 1/e of maximum), which indicates the com- 

(2) 

plete dephasing. As aforementioned, |xi;zz(w,t)| signi- 
fies the spin relaxation process. Thus we estimate that 
the spin relaxation time is about 10 fs, which is con- 
sistent with our previous results for time-resolved linear 



(magneto-)optics based on \x2y\ and \Xzz\- For the 
charge dynamics, we see a different scenario. In Fig. 

(2) 

1(b), \xzzz(u),t)\ is plotted as a function of time t. One 

(2) 

sees that \Xxzz\U, t) \ nearly spends a similar time to reach 
its maximum as \Xxzz{^>, t)\ does, but the subsequent de- 
cay occurs more sharply and strongly. Around 5 fs, the 

(2) 

value of |xzzz(w, t)\ is already close to the equilibrium 
data. This means that the dephasing already becomes 
strong for the charge dynamics before the spin dynam- 
ics. If one compares Fig. 1(a) with 1(b), one sees a 
clear difference between spin and charge dynamics. Ba- 
sically the spin dynamics lasts about twice as long as the 
charge dynamics. This has an important consequence 
as it demonstrates the spin memory effect: though the 
charge dynamics finishes, the spin dynamics is still alive, 
which is crucial for future applications. The main differ- 
ence of the time resolved nonlinear response compared 
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to the linear one, which is particularly evident for the 
magnetic dynamics, consists in an additional 'bunching' 
of the structures resulting from the simultaneous pres- 
ence of u> and 2uu resonances in Eqs. (2) and (3). 

In order to get a handle on the microscopic origin of 
the observed magnetic dynamics, we try to investigate 
the effect of the on-site exchange coupling J. We re- 
duce J to Jo/10. The corresponding time-dependences of 

(2) (2) 

\Xxzz(w,t)\ and IxLz^,^)! are shown in Figs. 1(c) and 

(2) 

1(d), respectively It can be seen that \Xxzz(u,t)\ first 

comes up within 2 fs. After that a recurrence appears 

with a rather large amplitude. Compared with Fig. 1(a), 
(2) 

\Xxzz{w, t) \ oscillates with a longer 'bunching' period and 
the loss of coherence is weaker. We estimate that the re- 
laxation time is about 10 fs but the period is nearly twice 
as long as that in Fig. 1(a). This demonstrates that the 
decrease of exchange interaction prolongs the period of 
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FIG. 1. The femtosecond time evolution of the nonlin- 
ear magneto-optical and optical responses, |x!tiz(w,t)| and £ 
|xlzz(w, i) I , respectively. Here the initial excited state is pre- 52- 00 
pared 2 eV above the ground state with a Gaussian broaden- Jj 
ing as wide as 20 eV. In (a) and (b) , a set of generic parameters 
of nickel is used; in (c) and (d) , the exchange interaction J is 
reduced to Jo/ 10 while the rest of the parameters is kept un- 
changed. The spin dynamics is delayed with respect to charge 
dynamics, which indicates the existence of a spin memory ef- 
fect. Additional 'bunching' occurs, which has not been found 

(2) 

for the linear pump-probe responses of Ni. \x x Jz(^,t)\ oscil- 
lates with a larger period. 



oscillations. For the charge dynamics, the change going 
from J to Jo/ 10 is relatively small. This can be seen 

(2) 

from Fig. 1(d) where we plot \xzzz(w,t)\ as a function of 
time t. Comparing Figs. 1(b) and 1(d), one finds that 
the overall variation of \xzzz{w,t)\ with time is nearly 
identical. This is understandable since the exchange in- 
teraction acts more directly on the spin degree of freedom 
by changing the spin dependence of the electronic many- 
body states microscopically. Consequently, the spin dy- 
namics will be affected more strongly than the charge 
dynamics. However electrons even with different spin- 
orientations play a similar role for the charge dynamics. 
Thus, the charge dynamics is basically independent of the 
spin state. That is why the exchange interaction does not 
affect the charge dynamics significantly. 

Next, we wish to gain some physical insights into the 
effects of the experimental constraints, such as the spec- 
tral width of the excited state distribution, on the spin 
dynamics. In Fig. 2, we perform a detailed compari- 
son between those experimental observables and intrinsic 
quantities (i.e., \S z (t)\ , |iV(t)| ). The initial excited state 
is prepared 2 eV above the ground state with a Gaussian 
broadening of now only 0.2 eV, which simulates a narrow 
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FIG. 2. Femtosecond time evolution of the intrinsic spin 
(a) and charge (b) dynamics, which is not directly observable, 
in comparison with the magneto-optical (c) and optical (d) re- 
sponse functions. |S z (t)| and |-/V(i)| , (a) and (b) show the 
clearly different behavior between spin and charge dynamics. 
Here the initial state is also prepared to be 2 eV above the 
ground state, but the spectral width is 0.2 eV. The other pa- 
rameters are taken as the generic parameters of Ni. The time 
interval is [0,40 fs] for (a) and (b). The calculated intrinsic 
((a), (b)) and experimentally accessible ((c), (d)) quantities 
show a strong dependence on the laser spectral width (corre- 
sponding to the initial excited state width). Comparing with 
Fig. 1, we notice that for the same set of generic parameters, 
the decrease of the laser spectral width yields a slower decay 
of the spin and charge response. 
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laser spectral width, in contrast to the previously used 
width of 20 eV. This allows us to see the effect of the 
laser spectral width clearly. All the parameters for the 
Hamiltonian are the generic set of parameters of Ni. In 
Figs. 2(a) and 2(b), we show the results up to 40 fs. The 
abscissa and ordinate in Figs. 2(a) and 2(b) denote the 
real and imaginary parts of the intrinsic spin and charge 
dynamics, S z (t) and N(t), respectively. Note that these 
quantities cannot directly be observed and are only theo- 
retically accessible. The arrows refer to the time direction 
and the centers are the final positions of S z (t) and N(t). 
One notices that the spin dynamics needs about six cycles 
to reach its final value while for the charge dynamics only 
three cycles are needed. This demonstrates again that 
the spin dynamics is delayed with respect to the charge 
dynamics. We find that comparing those intrinsic quanti- 
ties one can see a clear difference between spin and charge 
responses, which is not blurred by the details of the ex- 
perimental conditions. Nevertheless, from our results, 
one can still identify the differences also in the nonlinear 
optical and magneto-optical responses. We show the re- 
sults in Figs. 2(c) and 2(d). In Fig. 2(c), \Xxzz(u,t)\ is 
plotted as a function of time up to 80 fs. One may notice 
that a decay occurs within 30-40 fs, which is consistent 
with our previous linear results obtained from \xxy | and 
\Xzz\ f|- SHG probes more bands, which causes more 
'beating' and 'bunching'. This might then result in an 
effectively slower response than that seen in linear pump- 
probe experiments, in addition to the slowing down which 
results from the narrower bandwidth at interfaces. For 

( 2) 

|Xzzz(w,i)| we present the results in Fig. 2(d). Compar- 

(2) (2.) 

ing with \xxJz(u,t)\ in Fig. 2(c), \xzzz{u,t)\ drops even 
more sharply. One can see a clear drop at 18 fs, which 
sets its relaxation time for charge dynamics. A delay 
about 10 fs between spin and charge responses is found 
for this specific set of parameters. Comparing Figs. 1(a) 
and 2(c), one can immediately see that a narrow initial 
state distribution width prolongs the relaxation process. 
In Fig. 1, one knows that the relaxation process basi- 
cally finishes within 10 fs, but here the relaxation time is 
around 30-40 fs. 

Finally, as our previous studies already showed |(||7j for 
the linear pump-probe calculations, the band structure 
will also influence the relaxation process. Its effect is ac- 
tually very significant. Our results on | xi^ I an d IxS al- 
ready showed that the hopping integrals can modify the 
relaxation process strongly. Analogously, this will be re- 

(2) 

fleeted in the nonlinear optical responses \xxzz(to,t)\ and 

(2) 

\Xzzz(u, t)\ . In order to investigate the effect of the band 
structure, we reduce the hopping integrals to one-tenth 
of the original nickel hopping integrals while keeping the 
rest of parameters unchanged. Here, the initial excited 
state is also prepared 2 eV above the ground state with a 
Gaussian broadening of 20 eV. The results are shown in 
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FIG. 3. The femtosecond time evolution of the nonlin- 
ear magneto-optical and optical responses, \xxzz{u, t)\ and 
|xi««( w J *)l i respectively. The same initial condition has been 
used as in Fig. 1 except for the hopping integrals which are 
taken to be one-tenth of the original ones. It turns out that 
the change of band structure has a significant impact on spin 
and charge dynamics. It slows down the dynamics by modi- 
fying the envelopes of the response functions [x»*a(w,*)| and 

\xi 2 k",t)\ . 

(2) 

Fig. 3. One finds that the change of both \\ X zz{uJ^t) \ and 

(2) 

\Xzzz(u,t)\ with time is very different from the previous 

(2) 

ones. From Fig. 3(a), one notices that \xxzz(w,i)\ in- 
creases sharply within 5 fs. The strong oscillation lasts 
about 20 fs, where no clear decay can be seen. The out- 

(2) 

line of \xxzz (w, t) \ ranging from fs to 20 fs forms a broad 
peak. After 20 fs, dephasing occurs, but the envelope of 

(2) 

\Xxzz{oJ,t)\ decays only slowly. Comparing Figs. 1(a) 
and 1(b) with Figs. 3(a) and 3(b), respectively, one sees 
that the reduction of the hopping integrals slows down 
both the spin and charge dynamics considerably. The 

(2) 

envelope of \Xxzz{v,t)\ now decays on the time scale of 
30-40 fs. The appreciable fast oscillation in the long- 
time tail survives beyond 80 fs. Comparing Figs. 3(a) 
and 3(b) with Figs. 2(c) and 2(d), respectively, indi- 
cates that the elementary oscillations are more rapid. 
Thus from the comparison of all the three figures, we 
conclude that the reduction of the excited states dis- 
tribution width affects the elementary oscillation more 
strongly than the envelope of the spin and charge re- 
sponses while decreasing the hopping integrals does the 
opposite and mainly slows down the response envelope 

(2) (2) 

of \Xxzz(oJ,t)\ and |xWjz(w, £)| . It is remarkable that in 
all cases the spin memory effect persists. In addition, we 
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would like to point out that the change of band structure 
imposed by spin-orbit coupling as well affects the spin 
dynamics observed in nonlinear magneto-optical experi- 
ments. For linear time-resolved magneto-optical experi- 
ments, this has been demonstrated theoretically (see Ref. 
HQ). We found that spin-orbit coupling had to be raised 
to an unphysical large value of 1 eV (the generic value 
of spin-orbit coupling in nickel is 0.07 eV) in order to be 
relevant for the spin dynamics. Thus we do not expect 
any major influence from spin-orbit coupling on the 10 fs 
spin dynamics of transition metals in the valence band. 



IV. CONCLUSIONS 

In conclusion, we performed an exact-diagonalization 
scheme to study the spin dynamics in ferromagnetic 
nickel on the femtosecond time scale. This study is exclu- 
sively devoted to the nonlinear optical responses. At first, 
we checked that our Hamiltonian can give a correct ferro- 
magnetic ground state and a correct order of magnitude 
of the magnetic moment, which is larger than its value in 
bulk materials, as is physically expected. For a generic 
set of Ni parameters, the intrinsic speed limit of spin dy- 
namics is about 10 fs, which is on a similar time scale as 
revealed by \xxy\ although additional 'bunching' struc- 
ture appears for the nonlinear time-resolved magneto- 
optical response. Our theory clearly yields the memory 
effect of the spin dynamics also in nonlinear optics. It is 
found that the spin dynamics is delayed with respect to 
the charge dynamics. The spin dynamics survives even 
when the charge dynamics ceases. This is very impor- 
tant for future applications, such as ultrafast magnetic 
gates. We also examined the effect of exchange inter- 
action. It is found that a decrease of the exchange in- 
teraction prolongs the relaxation process. Especially we 
noticed that for a smaller exchange interaction J the spin 
response 'bunches' with longer periods. This again con- 
firms our earlier results that the observed spin dynamics 
results from the spin dependent ( J-dependent) dephas- 
ing in the excited many-body states. In that sense, SHG 
and NOLIMOKE reveal the same fundamental physics 
as time-resolved linear (magneto)-optics. In order to re- 
veal more insight into the spin dynamics, we also cal- 
culated some intrinsic quantities, from which we see a 
clearer difference between spin and charge dynamics. In 
addition, we studied the laser width effect on the spin 
dynamics. It shows that a narrow spectral width clearly 



slows down the dynamics. This is understandable as, for 
a smaller spectral width, the number of populated states 
is smaller. The dissipation becomes weaker and eventu- 
ally the spin dynamics is prolonged significantly. Finally 
we studied the effect of the band structure on the spin 
dynamics. We found that the reduction of the hopping 
integrals slows down both the spin and charge dynam- 
ics by modifying the envelope of the response functions 

(2) (2) 

\Xxzz(oJ,t)\ and |x*«a( w >*)l • This is different from the 
effect of a narrower distribution width, where one basi- 
cally prolongs the oscillation periods. It is important to 
note that materials with a small hopping integral corre- 
spond to nanostructure materials, clusters, quantum well 
states, magnetic insulators (such as oxides), defects, and 
impurities. Thus it is expected that with the presence 
of those nanostructure materials, the observed dynamics 
will be significantly slower. Naturally, this is very mean- 
ingful for the design of materials for ultrafast magnetic 
devices. Finally and mostly importantly, the processes 
which we investigated all exhibit the spin memory effect, 
a crucial property for future applications. 
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